Correlated fermions on a checkerboard lattice 
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A model of strongly correlated spinless fermions on a checkerboard lattice is mapped onto a quantum FPL 
model. We identify a large number of fluctuationless states specific to the fermionic case. We also show that 
for a class of fluctuating states, the fermionic sign problem can be gauged away. This claim is supported by 
numerical evaluation of the low-lying states. Furthermore, we analyze excitations at the Rokhsar-Kivelson 
point of this model using the relation to the height model and the single-mode approximation. 
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The interplay between quantum and geometric frustration 
can result in many unusual properties. From that point of 
view, strongly correlated spinless fermions on a checkerboard 
lattice present an interesting case. At certain fillings, fraction- 
ally charged excitations have been predicted for the case of 
large nearest-neighbor repulsion (U. The subject is not purely 
academic since the checkerboard lattice is a projection of a py- 
rochlore lattice onto a plane. There is experimental evidence 
that electrons in pyrochlore lattices can be strongly correlated 
J2]. We shall show that at half filling, this problem can be 
mapped onto a quantum fully packed loop (FPL) model - an 
analog of the quantum dimer model (QDM) - and discuss a 
number of consequences. QDMs, originally introduced in the 
context of quantum magnetism yD, became a focus of recent 
interest following the discovery of a gapped liquid phase on 
a triangular lattice U In particular, it has been established 
that a gapped phase with deconfined excitations exists in 2D 
for QDMs on non-bipartite lattices JH El l^l Q] while on bipar- 
tite lattices, such as a square lattice, systems typically crystal- 
lize into a phase with a broken translation/rotation symmetry. 
The liquid phase is "shrunk" into a quantum critical point with 
gapless excitations HI Ell ■ In both cases, an effective gauge 
theory is a U(l) theory for such a critical point and a Z2 theory 
for a deconfined phase 1 111 . Several related models, such 



as the quantum six/eight-vertex models 1121 Il3l Il4l 11511 have 



been shown to conform to the same dichotomy: a model with 
orientable loops is critical and is described by a U(l) gauge 
theory. 

The aforementioned models share one important feature: 
the matrix elements connecting various states of the low- 
energy Hilbert space are all non-negative. The Hilbert space 
separates into different sectors so that all states within a sector 
are connected by the quantum dynamics while states belong- 
ing to different sectors are not. Then, by Perron-Frobenius 
theorem, the ground state (GS) of each sector is unique and 
nodeless 13, El- Much less is known about models with non- 
Frobenius dynamics which results in some negative matrix el- 
ements. A striking difference between Frobenius vs. non- 
Frobenius QDMs on the kagome lattice was found in [0, U3l : 
while the conventional QDM exhibits a gapped Z2 topologi- 
cal phase, its counterpart with different signs has an extensive 



GS degeneracy. A different non-Frobenius model of lattice 
fermions ill 811 has also been found to have an extensive GS 
degeneracy. 

The starting Hamiltonian for spinless fermions is: 
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the summation is over neighboring sites of a checkerboard lat- 
tice. We assume that < t <C V - the limit of strong correla- 
tions. At half-filling which is studied here, the potential term 
is minimized whenever there are two fermions on each "planar 
tetrahedron", i.e., a crisscrossed square in Fig.[TJa) (the tetra- 
hedron rule, Il9ll ). Configurations satisfying this rule can be 
represented by dimers on a square lattice connecting the cen- 
ters of the planar tetrahedra, with the particles sitting in the 
middle of the dimers (Fig.QJa)) lf20h . Hence, the states satis- 
fying the tetrahedron rule are represented by non-intersecting 
FPL on the square lattice. Different loop configurations are 
orthogonal to each other as they correspond to different loca- 
tions of fermions (in the tight-binding limit, the overlap be- 
tween the corresponding wavefunctions is neglected). 

We now turn our attention to the kinetic term in Hamilto- 
nian (Q]i. It creates configurations in which the tetrahedron 
rule is violated. But since we are interested in the low-energy 
physics we may eliminate these configurations and consider 
instead an effective ring-exchange Hamiltonian 12111 that acts 
within the FPL Hilbert subspace. To lowest non-vanishing or- 
der in t/V, such effective Hamiltonian becomes: 



{m,g> 

{m,B> 

ee 9 J2 (\A)(A\ + \A)(A\ \B)(B\ \B)(B\ 
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where the sums are performed over all polygons of perimeter 
six and g = 12t 3 /V 2 . Particles are located in the middle 
of the dimers. The different signs result from the number of 
permuted fermions which is either even or odd. 
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FIG. 1: (a) - A configuration of fermions satisfying the "two per 
tetrahedron" rule and its representation in terms of FPLs. A small- 
est "fiippable plaquette" is shaded (see text for details), (b-d) - 
Fragments of possible "frozen" states: (b) is fiuctuationless under 
any local ring exchanges while (c) and (d) are fiuctuationless only 
for fermionic ring exchanges, (e) - The long-range ordered "squig- 
gle" state maximizing the number of three-dimer ring exchanges 
(shaded), (f) - Analog of a plaquette phase, with three dimers pre- 
dominantly fluctuating around shaded rectangles, shown dimers are 
not participating in these exchanges. 



Symmetries and Conserved Quantities - The first impor- 
tant point is that this model is a variation of the quantum six- 
vertex (6V) model Il2lll4ll5ll 029I1 , However, there are two 
important differences resulting from the fermionic statistics. 
Firstly, only ring exchanges involving odd numbers of parti- 
cles are allowed; the contributions from clock- and counter- 
clockwise ring moves of an even number of fermions cancel. 
In particular, the smallest resonating plaquette has perimeter 
six, rather than four, as captured by Eq. (JJJl. Secondly, the 
sign of the ring-hopping term depends on the occupancy of 
the middle bond. Both differences have important implica- 
tions: the first one - for the ergodicity of the quantum dy- 
namics, while the second one leads to a sign problem in the 
quantum Monte Carlo dynamics. 

Let us now discuss symmetries and conservation laws per- 
taining to our model. Firstly, since its Hilbert space is equiva- 
lent to that of a 6 V model, it has a height representation 122 ] . 
The quantum dynamics can also be described in the height 
language lE^lE^tl . The immediate consequence is that no such 
local dynamics can change the global slopes n x {y) and n y (x) 
in x and y directions, which therefore form a set of conserved 
numbers. Additionally, there is an accidental symmetry in our 
model, specific to the quantum dynamics (|2j: it conserves the 
total numbers of both vertical and horizontal dimers (more- 
over, it conserves the number of horizontal dimers in even 
and odd rows separately, and similarly for the vertical dimers). 
Hence, all states can be classified with respect to these quan- 
tum numbers which are different from the aforementioned 
conserved slopes. Notice that these additional conservation 
laws will not survive if higher-order ring-exchange terms are 
added to the Hamiltonian (0). 

Sign Problem - The sign problem manifested in the op- 
posite signs of the A «-> A and B <-> B terms in Eqs. (|2]i can 
be avoided in certain (but not all) cases. Notice that B <-» B 
processes do not change the loop topology (or their number), 
while the A <-> A flips always do, with the possibilities pre- 



sented schematically in Fig.[2|a). Let us consider the coarse- 
grained surface with an even number of fermions remaining 
fixed at the boundary ("fixed" boundary conditions (BC)). 
Then fermion configurations are represented by closed loops 
as well as arcs terminating at the boundary. With no fermions 
at the boundary, there are closed loops only (Fig.[2{b)), which 
we orient as follows: (i) color the areas separated by the loops 
white and grey, with white being the outmost color; (ii) orient 
all loops so that the white regions are always to the right. In 
the presence of arcs (Fig. |2jc)) we can choose how to close 
them on the outside without intersections (the outside region 
has no dynamics). Letting white be the color at infinity, we 
orient the loops as described above. We now notice that the 
relative signs resulting from the A <-> A flips are consistent 
with multiplying each loop configuration by i l (—i) r , where 
r (I) is the number of the (counter-) clockwise winding loop 
Hence, by simultaneously changing the sign of the A <-> A 
flips and transforming the loop states \C) — > i 1 ^' (— i) r ^ |£), 
we cure the sign problem thus making the system effectively 
bosonic. (In fact, this fix can be related to a more general class 
of gaugeable signs discussed in iflrll .) 




FIG. 2: (a) Three different actions of the effective Hamiltonian on 
the topology of loop configurations, (b, c) Representation of config- 
uration as fully-packed directed loops. 
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FIG. 3: GS energy and energies of lowest excited states of the effec- 
tive Hamiltonian in sectors with different slopes (Kx, 0). Left panel: 
Data from an exact diagonalization of Haf on a 72 site cluster where 
the Fermi sign is taken into account. Right panel: Same data from a 
calculation with excluded Fermi signs. 

Note that this construction need not work for the periodic 
BC: firstly, only even-winding sectors on a torus allow two- 
coloring, and secondly, even for even windings, it is possi- 
ble to dynamically reverse the coloring while returning to the 
same loop configuration, resulting in non-trivial state count- 
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ing. It appears though that for the periodic BC on even tori 
(preserving the bipartiteness of the lattice), the lowest-energy 
states belong to the sector where such transformation works, 
as confirmed by the exact diagonalization results (Fig. [3]). The 
presented non-local loop-orienting construction is restricted 
to the ring-exchange processes of length six. 

Rokhsar-Kivelson (RK) Point - It has been shown that at 
half-filling the Hamiltonian (f2l) has an ordered GS and frac- 
tionally charged, weakly confined excitations B25I1 . In the 
spirit of 13], we extend our model's parameter space in order 
to study its phase diagram and consider the Hamiltonian 

Kn = E [-9 i\ A )&\ + \A)(A\ \B)(B\ \B)(B\) 



GSs are not related to any height constraints - in fact, the 
states of type shown in (Fig.[Tfd)) include both the maximally 
tilted state (Fig.QIb)) and the maximally flat state (Fig.[TJc)). 
"Freezing" of the latter state is especially interesting since it is 
supposed to dominate the ordered phase of the bosonic model 
for/x/g < —0.374 J 1311 by maximizing the number of the flip- 
pable plaquettes 13 ill . This observation puts into doubt the 
existence of such a phase in the fermionic model. 

In addition to the frozen states, there are also liquid-like 
GSs, as was the case in the original RK construction. These 
states can be found exactly for the dynamics given by Eqs. (f2|- 
@]i with the help of the aforementioned sign-fixing transfor- 
mation. (While this transformation changes the sign of the 
A <-» A flips, it does not affect the diagonal terms in Eqs. ([3]).) 



-H(\A){A\ + \A)(A\ + \B)(B\ + \B){B\)] (3) Hence, for "fixed" BC, the state |0) oc£ {£} 



where we use the short-hand notations introduced in Eq. (0). 
At the RK point g = fx, this Hamiltonian can be written as a 
sum of projectors: 



gJ2[{\A) + \A))((A\ + {A\) 
p 

+ (\B)-\B))((B\-(B\)}. (4) 

A (non-unique) GS of Hamiltonian (O is a state that is anni- 
hilated simultaneously by all projectors. There are two kinds 
of such states: frozen states and liquid states. It turns out that 
the fermionic nature has strong consequences for both kinds. 
All bosonic frozen states identified in 111311 remain frozen since 
they are unaffected by any local ring exchange. In the height 
language, these states maximize the slope in either x or y di- 
rection, making any local dynamics impossible. Fixing this 
slope in an L x L system, one still has 2 L possible arrange- 
ments of steps in the other direction which results in 4 x 2 L 
such states [30]. However, there are also new additional GSs 
(see Fig. QI c id)) which remain unaffected by any fermionic 
dynamics but would not be so in the bosonic case. This is 
due to the fact that only odd-fermion ring exchanges are al- 
lowed here. For the types of states shown in Fig. Q2c,d), this 
is easy to see: any ring exchange is an alternating sequence 
of occupied and vacant bonds. The subsequence of vacant 
bonds maps onto a ID walk between the lines and must con- 
tain an even number of steps in order to close. The states 
of this type are obtained by translating a single line like in 
Fig- Old), e.g., a line starting at the left bottom corner (0, 0) 
and reaching the top boundary at (x, L) while going only up 
and to the right (we assume free BC) can do it in ( L x x ) dis- 
tinct ways leading to 4 x J^t=o ^f) — * 4 L+1 ( for lar 8 e 
L) such frozen fermionic states. There is also a number of 
other frozen fermionic configurations that we could pinpoint 
(but not systematically count). Hence it remains an open is- 
sue whether our RK point is characterized by an extensive 



GS entropy and whether such entropy observed in IU7L 11811 



can be explained in this manner. However, it is very sugges- 
tive that the explanation lies in a more restrictive nature of 
the fermionic dynamics. We stress that these "new" frozen 



is automatically annihilated by all projector terms in H e ff and 
therefore is a GS. Here the sum is taken over all loop config- 
urations £ that are connected to each other by the quantum 
dynamics. Note that there are many disconnected sectors. 

Excitations at the RK point - The quantum dynamics of 
the Hamiltonian H e s and similarly H' e{{ can be described in 
terms of a height model 123l 12411 . The associated conserved 
quantities k x and K y , as well as their insensitivity to a constant 
shift of the height field imply gapless hydrodynamic modes 
with uj(k) ~ k 2 as k — ► 0. The liquid state of the FPL model 
at the RK point corresponds to the rough phase of the height 
model, in which case the modes can be identified as the so- 
called resonons at wavevecto r (tt , tt) Oland the equivalent of 
the piOns, here at Q = (0, 0) l^lH lm 

For constructing gapless excitations in our model, we adopt 
the single-mode approximation which has been successfully 
used for the QDM on the square lattice fl Let |0) be 
the aforementioned "liquid" GS at the RK point. The opera- 
tor (r) creates (annihilates) a dimer at position r in direc- 
tion f. The density operator n?(r) = df (r)dj (r) has the 
Fourier transform rif(q) = ^ r n+ (r) exp(iq ■ r). We use 
the operators rif (q) to construct the states |q, f) = (q)|0) 
which are for q ^ orthogonal to |0). The excitation 
energies are bounded by E(q, f) < /(q)/SV r (q), where 
/(q) = (0|[nf(-q), [F^ ff ,nf(q)]]|0) is the the oscillator 
strength and S TT (q) = (0|nf (— q)n,f (q)|0) is the structure 
factor. In order to calculate /(q), we observe that the den- 
sity operators commute with the potential energy term, hence 
only the kinetic energy term contributes. By using the com- 
mutation relation [df,nf] = =fdf repeatedly, we compute 
/(q) due to the resonating terms for f = x. The upshot is 
that to leading order /(k) ~ (k x t) 2 (similar type of ex- 
pression as the one for the RK model yfl) with k = q — Q at 
Q = (0, 0), (tt, tt), (0, tt). The second ingredient, the struc- 
ture factor S xx (q) is the Fourier transform of the dimer den- 
sity correlation function. We use the expression given in 
Ref. II27I1 for it. Interestingly, S xx (q) ^ except along the 
direction q = (q x , tt) where it vanishes with the exception 
of (tt, tt). At (0, tt) both /(q) and S xx (q) are zero, but their 
ratio remains finite. In addition S xx (q) shows no singulari- 
ties. The FPL model differs from the hardcore dimer model 
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(describing our model at quarter-filling) for which S xx (q) di- 
verges logarithmically at Q = (it, 0) 12811 . The difference is 
due to a slower algebraic decrease with distance of the cor- 
relation function for hardcore dimer covering. We have also 
verified the above result for our structure factor S X x(<l) by 
means of Monte Carlo simulations. 

We conclude with several observations and open questions 
addressing the possible phase diagram of the model defined by 
Eq. y). The bosonic version studied in [13] has three phases 
depending on the parameter fj,/g: the two "fiat" phases, i.e., 
the Neel phase 1 1311 and the plaquette phase (also found in 
lfl4l , I27I1 ) as well as a maximally tilted frozen phase, with the 
RK point perched between the latter two phases. We have al- 
ready argued that some of the flat states are actually frozen for 
fermions and hence cannot be used to define the fluctuation- 
stabilized phases for /i/g < 1. The nature, and number of 
such phases are an open question, but let us attempt an anal- 
ysis. Firstly, the maximally-flat (Neel, DDW) phase shown 
in Fig. [TJc) appears to be replaced by the somewhat analo- 
gous "squiggle" phase (Fig. [TTe)) ll25ll . In the bosonic case, 
the region of —0.374 < fi/g < 1 corresponded to the plaque- 
tte phase, while its direct fermionic analog does not appear 
to be present anywhere in the fermionic phase diagram. (No- 
tice, however, that unlike the Neel and the plaquette phases, 
both states shown in Fig. [TJ e T) break both translation and 
rotation symmetries simultaneously.) We rule this phase out 
based on the fact that the squiggle phase breaks the symmetry 
between the numbers of the vertical and the horizontal dimers 
(which are separately conserved by the quantum dynamics), 
while the "fermionic plaquette" phase does not. The results 
of the exact diagonalization on the small samples indicate that 
this symmetry remains broken all the way up to fi/g = 1 - 
the RK point, with the ratio of 2/3 being consistent with the 
"squiggle" phase. This observation is also strongly disfavor- 
ing a critical liquid phase to the left of the RK point, hence 
the fermionic RK point in this model is likely to be an isolated 
quantum critical point just as it is for the bosonic model. The 
gapless modes identified above are in agreement with such a 
scenario. The details of the phase diagram and its sensitivity 
to the higher order ring exchanges is a subject for future work. 
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The mapping between FPL and 6V is done by subdividing the 
square lattice into even and odd sublattices and orienting occu- 
pied bonds from odd to even sites, with the opposite orientation 
for vacant bonds. 

This is unlike the QDM case with only 4 such frozen, or stag- 
gered states: maximizing the slope in one direction there auto- 
matically fixes it in the other. 

This is the conventional anti-ferroelectric phase of the six- 
vertex model referred to as the Neel phase in Jl3ll and as DDW 
phase in 0GJ] 



